Subsystem Codes 



Salah A. Aly, Andreas Klappenecker and Pradeep Kiran Sarvepalli 



Abstract — We investigate various aspects of operator quan- 
tum error-correcting codes or, as we prefer to call them, 
subsystem codes. We give various methods to derive subsystem 
codes from classical codes. We give a proof for the existence 
of subsystem codes using a counting argument similar to the 
quantum Gilbert- Varshamov bound. We derive linear program- 
ming bounds and other upper bounds. We answer the question 
whether or not there exist [[n, n — 2d + 2,r > 0, d]] q subsystem 
codes. Finally, we compare stabilizer and subsystem codes with 
respect to the required number of syndrome qudits. 

I. INTRODUCTION 

Quantum error-correcting codes are seen as being in- 
dispensable for building a quantum computer. There are 
three predominant approaches to quantum error-correction: 
stabilizer codes, noiseless systems, and decoherence free 
subspaces. Recent advances in the theory of quantum error- 
correction have shown that all these apparently disparate 
approaches are actually the same. This unification goes under 
the name of operator quantum error-correction codes [1], 
[7]-[ll], though we will prefer to use the shorter and more 
descriptive term subsystem codes. Subsystem codes provide 
a common platform for comparing the various different types 
of quantum codes and make it possible to treat active and 
passive quantum error-correction within the same framework. 
Apart from the fact that subsystem codes give us more 
control over the degree of passive error-correction, there have 
been claims that subsystem codes can make quantum error- 
correction more robust and practical. For example, it has 
been claimed that subsystem codes make it possible to derive 
simpler error recovery schemes in comparison to stabilizer 
codes. Furthermore, it is conjectured that certain subsystem 
codes are self-correcting [1]. 

Subsystem codes are relatively new and promise to be a 
fruitful area for quantum error-correction. Until now, there 
are few concrete examples of such codes and even fewer 
systematic code constructions. Little is known about the 
parameters of subsystem codes, so it is difficult to judge 
the performance of such codes. 

In a recent work [6], we derived a character-theoretic 
framework for the construction of subsystem codes. We were 
able to show that there exists a correspondence between the 
subsystem codes over F q and classical additive codes over 
F q and ¥ q 2 . In this paper, we investigate basic properties of 
subsystem codes, establish further connections to classical 
codes, and derive bounds on the parameters of subsystem 
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codes. We report first results on a fair comparison between 
stabilizer codes and subsystem codes. 

The paper is structured as follows. After a brief introduc- 
tion to subsystem codes in Section||l] we recall some results 
about the relations between subsystem codes and classical 
codes. Then we give some simple constructions of subsystem 
codes which parallel the common constructions of stabilizer 
codes. In Section!!!]] we give a nonconstructive proof of the 
existence of subsystem codes. In Section lTVl we derive linear 
programming upper bounds on the parameters of subsystem 
codes. For pure subsystem codes (to be defined later) we 
can also derive analytical upper bounds which resemble 
the quantum Singleton and Hamming bounds. Armed with 
these results on bounds we make a rigorous comparison of 
stabilizer codes and subsystem codes, that makes precise 
when subsystem codes can do better. 

Notation: We assume that q is the power of a prime p 
and F q denotes a finite field with q elements. By qudit we 
mean a q-ary quantum bit. The symplectic weight of an 
element w = (xi, . . . , x n , y\, . . . , y n ) in F q n is defined as 
swt(w) = \{(xi,yi) ^ (0,0) | 1 < i < n}\. The trace- 
symplectic product of two elements u = (a\b),v = (a'\b') 
in Fq U is defined as (u\v) s = ti q / p (a' ■ b — a ■ b') , where x ■ y 
is the usual euclidean inner product. The trace-symplectic 
dual of a code C C F 2 q n is defined as = {v € F q n | 
(v\w) s = for all w G C}. For vectors x,y in F™ 2 , we 
define the hermitian inner product (x\y)h = XaLi x \Vi an d 
the hermitian dual of C C F£ 2 as C ±h = {x € F™ 2 | 
(x\y)h = for all y 6 C}. The trace alternating form of two 
vectors u 7 w in F™ 2 is defined as {u\v) a = tr q / p [((u\v)h — 
{v\u) h )/(f3 2 -(3 2q )}, where {f3,(3 q } is a normal basis ofF g2 
over Fq. If C C F™ 2 , then the trace alternating dual of C is 
defined as C Xa = {x g F™ 2 | {x\y) a = for all y € C}. 

II. SUBSYSTEM CODES AND CLASSICAL CODES 

Let H = C q ® C q • • ■ ® C q = C q ". An orthonormal 
basis for C q is B = {\x) \ x e F™}. The vector \x) = 
\xi) ® \x2) <8> • ■ • (8 \x n )- The elements of H are of the form 

v = v x \x) where v x £ C and \ v x\ 2 = L 

We define the following unitary operators on C 9 

X a \x) = \x + a) and Z b \x) = w tr ^ (fKE) |x), 

where u = e j27r / p . The set of errors £ = {X a Z b \ a,b € F q } 
form a basis for errors on a single qudit. Every error on a 
single qudit can be expressed as linear combination of the 
elements in £ . If we assume that the errors are independent 



on each qudit, we need only consider the error group E = 
{uo c e\ ® e 2 ® • • • <8> e n \ c G F p , a G £}, where each of the 
is a single qudit error. The weight of an error is the number 
of qudits that are in error. For further details on the error 
model and the actual structure of the error group we refer 
the reader to [5]. 

A quantum error-correcting code Q is a subspace in TL = 
C 9 " such that TL = Q © Q 1 -, where Q 1 - is the orthogonal 
complement of Q. In a subsystem code, the subspace Q 
further decomposes into a tensor product of two vector space 
A and B, that is, 

Q = A® B. 

The vectors spaces A and B are respectively called the sub- 
system and the co-subsystem of the code Q. The information 
to be protected is stored in the subsystem A, whence the 
name subsystem code. 

If dim A = K, dim£? = R and Q is able to detect all 
errors in E of weight less than d on subsystem A, then we 
say that Q is an {{n,K,R,d)) q subsystem code. We call 
d the minimum distance of the subsystem A or, by slight 
abuse of language, the minimum distance of the subsystem 
code Q (when the tensor decomposition Q = A ® B is 
understood from the context). We write [[n, k, r, d]] q for an 
((n, q k ,q r ,d)) q subsystem code. Sometimes we will say that 
an [[n, k, r, d]] q subsystem code has r virtual gauge qudits, 
which is simply another way of saying that the dimension of 
the co-subsystem is q r ; it should be stressed that the gauge 
qudits typically do not correspond to physical qudits. 

A. Subsystem Codes From Classical Codes 

We recall the following results from [6] which relate 
quantum subsystem codes to classical codes. 

Theorem 1: Let X be a classical additive subcode of ¥ 2n 
such that X ^ {0} and let Y denote its subcode Y = X n 
X ±a . If x = \X\ and y = \Y\, then there exists subsystem 
code C = A® B such that 

i) dim ,4 = q n /{xy) 1 / 2 , 

ii) dimB = {x/y) 1 / 2 . 

The minimum distance of subsystem A is given by d = 
swt{{X + X^°)-X) = swt(y ±s -X). Thus, the subsystem 
A can detect all errors in E of weight less than d, and can 
correct all errors in E of weight < \_(d — 1 ) / 2 J . 

Proof: See [6, Theorem 5]. ■ 

Remark 1: Recall that \X^*\ = q 2n /\X\. Therefore, the 
dimension of the subsystem A can also be calculated as 

dim^i = (\x x '\/\y\)V 2 . 

It is also possible to construct subsystem codes via codes 
over ¥ q 2 using ( • | ■ ) a , the trace alternating form [5] which 
gives us the following theorem. The proof can be found in 
[6, Theorem 6]. 

Theorem 2: Let X be a classical additive subcode of F™ 2 
such that X i {0} and let Y denote its subcode Y = X n 
X a . If x = \X\ and y = \Y\, then there exists subsystem 
code C = A® B such that 

i) dimA = q n /(xy) 1 / 2 , 

ii) dimB = {x/y) 1 / 2 . 



The minimum distance of subsystem A is given by 

d = wt((X + X ±a ) - X) = wt(F- L ° - X), 

where wt denotes the Hamming weight. Thus, the subsystem 
A can detect all errors in E of Hamming weight less than 
d, and can correct all errors in E of Hamming weight [(d — 
1)/2J or less. 

Proof: This follows from Theorem H an d the fact that 
there exists a weight-preserving isometric isomorphism from 
{¥ 2n ,(-\-) s ) and (F» 2 ,(-|-}a), see [5]. ■ 

Theorem|2]has the advantage that the weights of the codes 
over ¥ q 2 is measured using the usual Hamming distance. 

We are now going to derive some particularly important 
special cases of the above two theorems as a consequence. 
Before stating these results, we recall the following simple 
fact. 

Lemma 3: Let C\ and C% be two ¥ q -linear codes of length 
n. The product code C± x C2 — {(a\b) \ a S C\,b £ C2] 
has length In and its trace-symplectic dual is given by 

(Ci x C a ) x ' = Ci x Ci. 
Proof: If (a\b) eC 1 xC 2 and (a'\b') G x C u then 
tr 9/p (6 • a' - V ■ a) = 0; hence, x C {C\ x C 2 ) ±s . 
Comparing dimensions shows that equality must hold. ■ 

The first consequence uses the euclidean inner product, 
that is, the usual dot inner product on F™ to construct 
subsystem codes. In the special case of stabilizer codes, this 
yields the well-known CSS construction (for instance, see 
[3, Theorem 9]). 

Corollary 4 (Euclidean Construction): Let C,; C F™, be 
[n,ki] q linear codes where i G {1,2}. Then there exists an 
[[n, k, r, d]] q subsystem code with 

. k = n- (ki+k 2 + k')/2, 

. r = (fci + fc 2 - k')/2, and 

. d = min{ wt ( ( d n C 2 ) - 1 \ d ) , wt ( {Ci n Ci ) x \ C 2 ) } , 
where k' = dim Fg (Ci n C^) x (Cf n C 2 ). 

Proof: Let C = C x x C 2 , then by Lemma C x > = 
C£- x C^, and D = c n c^ = {d n C^) x (C 2 n ci). 
Again by Lemma|3] = (C 2 nC i L )- L x (C x HC^)- 1 . Let 
dim F? L> = k'. Then |C||D| = q ^+^+ k> and |C|/|D| = 
qk 1 +k 2 -k gy Theorem n the code C defines an [[n,n — 
(fci + fc 2 + k')/2, (fci + fc 2 — k')/2, d]] q subsystem code. The 
distance of the code is given by 

d = swt(L>^ \ C) 

= swt((Ca n ci) 1 - x (Ci n Ci) 1 - \ (Ci x C2)). 

The latter expression can be simplified to 

d = min{wt((C 2 n Ci) 1 - \ Ci),wt((Ci D C^) x \ C 2 )}, 

which proves the claim. ■ 
Setting C 2 = Ci in the previous construction simplifies 
the computation of the code parameters. Then we have an 

[[n,n - k - k' ,k - k' ,wt{{Ci nCi)^ \ Ci)}} q code, where 
k' = dimp 5 CiflCi. Therefore, any family of classical codes 
where the dimension of C\ n Ci and the minimum distance 
of the dual of C\ n Ci is known, will provide us with a 
family of subsystem codes. The codes that arise when C\ = 



C2 will also arise as a special case of the next construction. 



Corollary 5 (Hermitian Construction): Let C C F™ 2 be 
an F ? 2 -linear 



dime 



dimension k' 

[In, n — k — k ,k 



code such that D = C H C ±h is of 
D. Then there exists an 



fe',wt(D^\C)]] ? 



subsystem code. 

Proo/:- If C* is linear, then C ±a = C ±h by [6, 

and 



i—k —k 



L emma 18 ], It follows that q n /^\D\\C\ = q 
y/\C\/\D\ = q k - k '. Let d = wt{D ±h \ C). Then, by 
Theorem there exists an [[n,n — k — k',k — k',d]] q 
subsystem code. ■ 
The subsystem codes can be easily constructed with the help 
of a computer algebra system. The following example gives 
some subsystem codes constructed using MAGMA [2]. 

Example 1 (BCH Subsystem Codes): The binary subsys- 
tem codes in Table U were derived from BCH codes over F4 
via Corollary |3 

TABLE I 
BCH Subsystem Codes 



Subsystem Code 


Parent 
BCH Code 


Designed 
distance 


[[15, 1, 2, 5]] 2 
[[15,5,2,3]] 2 


[15,8, 6] 2 2 
[15,6, 7] 2 2 


6 
7 


[[17,8,1,4JJ 2 


[17,5,9j 2 2 


4 


[[21,6,3,3JJ 2 
[[21,7,2,3]] 2 


[21,9,7jj 2 2 
[21,8, 9] 2 2 


6 

8 


[[31,10, 1,511a 
[[31,20, 1,3]] 2 


[31,11, 11J 2 2 
[31,6,15] 22 


8 

12 



Codes constructed with the help of Corollaries 0] and 
13 will lead to F g -linear and F ? 2 -linear subsystem codes 
respectively. Though in some cases Corollary |4] can lead to 
F g 2-linear codes. So when we refer to a subsystem code as 
being F g -linear, it could be also ¥ q 2 -linear. In this paper, 
we will call a subsystem code that can be constructed 
with the help of Theorems \l\ and [2] and their corollaries, 
a Clifford subsystem code. Ten years ago, Knill suggested 
a generalization of stabilizer codes that became known as 
Clifford codes (because their construction uses a part of rep- 
resentation theory known as Clifford theory). Recently, we 
realized that a special case of Knill's construction leads to a 
very natural construction of subsystem codes. Clifford theory 
is the natural tool in the construction of these subsystem 
codes, whence the name. 

III. LOWER BOUNDS ON SUBSYSTEM CODES 

In this section we give a simple nonconstructive proof for 
the existence of subsystem codes. The proof is based on a 
counting argument similar to the quantum Gilbert- Varshamov 
bound for stabilizer codes [5]. We will need the following 
simple fact. 

Lemma 6: Let F g be a finite field of characteristic p. Let 
r and s be nonnegative integers such that p r+2s < q 2n . Then 
there exists an additive subcode X of ¥^ n such that \X\ = 



Proof: Let m denote the integer such that q = p m . 
We may regard F^™ as an 2nm-dimensional vector space 
over F p . Then (• | -} s is a nondegenerate skew-symmetric 
bilinear form on this vector space. Therefore, there exists a 
direct sum decomposition of F^™ = ¥ 2nm = Vi © ■ • • ffi V nm , 
where V k is a 2-dimensional subspace with basis {x k ,Zk} 
such that (xk \ xe) s = = (zk \ zi) s for 1 < k,£ < nm, 
(x k I z k ) s ^ 0, and (x k | z e ) s = if k ^ I. Then X = 
(zi, . . . , z r , av+i, z r -\-\, . . . , x r + s , z r + s ) is a code with the 
desired properties. ■ 
Theorem 7: Let ¥ q be a finite field of characteristic p. If 
K and R are powers of p such that 1 < KR < q n and d is 
a positive integer such that 

d-l , s 

E f tf- l) j {q n KR - q n R/K) < (p - l)(q 2n - 1) 

holds, then an ((n, K, R, > d)) q subsystem code exists. 

Proof: By Lemma|6] there exists an additive subcode X 
of F^" such that x = \X\ = q n R/K and y = \X n X x ' \ = 
q n /(KR); the resulting subsystem code has a subsystem 
of dimension q n /(xy) 1 ^ 2 = K and a co-subsystem of 
dimension (x/y) 1 ^ 2 — R. Therefore, the multiset X given 

by 



X 



(X + X ±S ) - X 



X is an additive subcode of 
¥ 2 q n such that \X\ = q n R/K 
and \XnX ±s \ = q n /(KR) 



is not empty. 

Thus, an element of X corresponds to a subsystem code 
C = A® B with dimA = K and dimB = R. The set 
difference (X + X ±s ) — X contains only nonzero vectors 
of ¥ 2n . We claim that all nonzero vector in F^ n appear in 
the same number of sets in X. Indeed, the symplectic group 
Sp(2n, F q ) acts transitively on the set F^" \ {0}, see [4, 
Proposition 3.2], which means that for any nonzero vectors 
u and v in ¥ 2n there exists r G Sp(2n, ¥ q ) such that v = tu. 
Therefore, u is contained in (X + X ±s ) — X if and only if 
v is contained in the element (tX + (tJ) 1s ) — tX of X. 

Since \(X + X ±s ) - X\ = q n KR - q n R/K, we 
can conclude that any nonzero vector of F^' 1 occurs in 
\X\(q n KR-q n R/K)/(q 2n - 
a nonzero vector and its F* 



T) elements of X. Furthermore, 
-multiples are contained in the 
exact same sets of X. Therefore, if we delete all sets from 
X that contain a nonzero vector with symplectic weight less 
than d, then we remove at most 

ElZlOtf-lY lvi (q»KR-q»R/K) 



p-1 



X\ 



- 1 



P 



r+2s and IX n X ± ° 



P 



sets from X. By assumption, this number is less than \X\; 
hence, there exists an ((n, K, R, > d)) q subsystem code. ■ 
The lower bound has important implications for comparing 
stabilizer codes with subsystem codes as we shall see in 
Section IVII Further, we obtain the following lower bound 
for stabilizer codes as a simple corollary, when R = 1 (see 
also [5]). 



Corollary 8 (GV Bound for Stabilizer Codes): Let F g be 
a finite field of characteristic p and 1 < K < q n a power of 
p. If 



d-l 



iy < (p - 1) 



Or" - 1) 



(q n K -q n /K) 



holds, then an ((n,K, > d)) q stabilizer code exists. 

A stronger result showing the existing of linear stabilizer 

codes was shown in [5, Lemma 31]. 

IV. UPPER BOUNDS FOR SUBSYSTEM CODES 

We want to investigate some limitations on subsystem 
codes that can be constructed with the help of Theorem ^ 
(or, equivalently, Theorem^. To that end, we will investigate 
some upper bounds on the parameters of subsystem codes. 

A. Linear Programming Bounds 

Theorem 9: If an ((n, K, R,d)) q Clifford subsystem code 
with K > 1 exists, then there exists a solution to the 
optimization problem: maximize 2~2jZi Ai subject to the 
constraints 

1) A = B a = 1 and < Bj < Aj for all 1 < j < n; 

n n 

2) a j = v nR / K > B i = tl^R; 



3) A, 3 = — - Kj(r)A r holds for all j in the range 

q n R ' — ! 

r— 

< j < n; 

KR ™ 

4) Bj- S = — — Kj(r)B r holds for all j in the range 

< j < n; 

5) Aj = Bf' for all j in < j < d and Aj < Bf 3 for 
all d < j < n; 

6) B-j = Aj 3 for all j in < j < d and Bj < Af" for 
all d < j < n; 

7) (p - 1) divides Aj, Bj, Af 3 , and Bf" for all j in the 
range 1 < j < n; 

where the coefficients Aj and Bj assume only integer values, 
and Kj (r) denotes the Krawtchouk polynomial 



^•(r) = £(-l) s (g 2 -iy- , 

\ s I V 7 — s . 

s=0 \ / \J / 

Proof: If an ((n, K, R, d)) q subsystem code exists, then 
the weight distribution Aj of the associated additive code X 
and the weight distribution Bj of its subcode Y = X f] 
X ±s obviously satisfy 1). By Theorem ^ we have K = 
q n /y / \X\\Y\ and R = ^\X\/\Y\, which implies \X\ = 
Y,Aj = q n R/K and \Y\ = j^Bj ; = q n /KR, proving 2). 
Conditions 3) and 4) follow from the Mac Williams relation 
for symplectic weight distribution, see [5, Theorem 23]. As 
X is an F p -linear code, for each nonzero codeword c in X, 
ac is again in X for all a in F p ; thus, condition 7) must hold. 
Since the quantum code has minimum distance d, all vectors 
of symplectic weight less than d in Y ±3 must be in X, since 
Y ±s —X has minimum distance d; this implies 5). Similarly, 
all vectors in X ±3 C X + X s of symplectic weight less 



than d must be contained in X, since (X + X s ) — X has 
minimum distance d; this implies 6). ■ 

We can use the previous theorem to derive bounds on the 
dimension of the co-subsystem. If the optimization problem 
is not solvable, then we can immediately conclude that a 
code with the corresponding parameter settings cannot exist. 

Perhaps one of the most striking features of subsystem 
codes is the potential reduction of syndrome measurements. 
Recall that an F g -linear [[n, fc, d]] q stabilizer code requires 
n — k syndrome measurements. On the other hand, an F g - 
linear [[n, k, r, d]] q Clifford subsystem code requires just n — 
k — r syndrome measurements. 

Poulin [11] asked whether we can have [[5, l,r > 0, 3]]2 
Clifford subsystem code. Of course, such a code would 
be preferable over the [[5, 1,3]] 2 stabilizer code. After an 
exhaustive computer search, he concluded that such a sub- 
system code does not exist. This result can be obtained 
very easily with the linear programming bounds. In fact, 
our investigations for small lengths revealed that not only 
a [[5, 1, r > 0, 3]] 2 code does not exist, but neither does any 
code with parameters given in the next example. 

Example 2: Theorem [9] shows that it is not possible to 
construct subsystem codes with r > and parameters shown 
in Table HD 

TABLE II 



Field 


Codes 


F 2 


[[4,2,r,2]] 2 , [[5,l,r,3]] 2 


F 3 


[[4,2,r,2]] 3 , [[5,l,r,3]] 3 , 
[[9,3,r,4]] 3 , [[9,5,r,3]] 3 , 
[[10, 6, r, 3]] 3 


F 4 


[[4,2,r,2JJ 4 , [[5,l,r,3jj 4 , 
[[9,3,r,4]] 4 , [[9,5,r,3]] 4 , 
[[10, 6, r, 3]] 4 



The previous example is motivated by the fact that one 
can improve upon Shor's [[9, 1,3]]2 quantum stabilizer code 
by allowing three additional gauge qubits, that is, there 
exists a [[9, 1, 3, 3]]2 subsystem code, see [11]. The practical 
relevance is that the 9 — 1 = 8 syndrome measurements that 
are required for Shor's code are reduced to 9 — 1 — 3 = 5 
syndrome measurements in the subsystem code. 

Since we allow nonbinary alphabets in this paper, a natural 
generalization of Poulin's question is whether one can find 
an [[n,n — 2d+ 2,r, d]] q subsystem code with r > 0. The 
above example shows that such subsystem codes with such 
parameters do not exist for certain small lengths and small 
alphabet sizes. 

We will fully answer this question in the subsequent 
sections. In the search for an answer to this problem, we 
were prompted to define the notion of pure subsystem codes. 
The notion of purity proved to be fruitful in deducing this 
and other results. 

B. Pure Subsystem Codes 

Let X be an additive subcode of F^' 1 and Y = Xf)X ±3 . 
By Theorem^ we can obtain an ((n,K,R,d)) q subsystem 



code Q from X that has minimum distance d = swt(y ±s — 
X). The set difference involved in the definition of the 
minimum distance makes it harder to compute the minimum 
distance. Therefore, we introduce pure codes that are easier 
to analyze. 

We say that the subsystem code Q is pure to d! if d' < 
swt(X). The code is exactly pure to d! if it is pure to d! 
but not to d! + 1; then swt(X) = d! . Any subsystem code 
is always exactly pure to d' = swt(X). We call Q a pure 
subsystem code if it is pure to d 1 > d; otherwise, we call 
Q an impure subsystem code. Pure codes do not require 
us to compute the minimum distance of the difference set 
Y ±s — X. We can compute the distance of the code as d = 
swt(F is ), which is comparatively simpler task though it is 
also computationally hard. 

The purity of codes over ¥ q 2 is defined in a similar way. 

Example 3 (Reed-Solomon Subsystem Codes): The non- 
binary subsystem codes given in Table |ni| are all pure 
and were derived from primitive narrowsense Reed-Solomon 
codes over ¥ q 2. It is curious that the distance of many of 



TABLE III 
Reed-Solomon Subsystem codes 



Subsystem Codes 


Parent 
RS Code 


[[15, 1, 10,3|| 4 
[[15,1,8,3]] 4 
[[15,1,6,3]] 4 
[[15,2,5,3]] 4 


[15, 12, 4J 4 2 
[15,11,5] 4 2 
[15,10,6] 4 2 
[15,9,7] 4 2 


[[24,1,17,411s 
[[24,2,10,4]] 5 
[[24,4, 10,4]] 5 
[[24, 16,2,4]] 5 
[[24,17,1,4]] 5 
[[24,19,1,311s 


[ 


24, 20, 5 
24, 16, 9 
24, 15, 1C 
24,5,20 
24,4,21 
24,3,22 


5 2 
5 2 
] 5 2 
5 2 
5 2 
5 2 


[[48,1,37,6]] 7 
[[48, 2, 26, 6]] 7 


[ 


48,42,7 
48,36, IS 


7 2 

It 2 



these subsystem codes is equal to q—1. We conjecture that, 
in general, the distance of a subsytem code constructed from 

2 _ 

a Reed-Solomon code over ¥ q q2 cannot exceed q — 1. 
C. Upper Bounds for Pure Subsystem Codes 

In this subsection, we establish a number of basic results 
concerning pure subsystem codes. The next lemma is a 
key result that associates to a pure subsystem code a pure 
stabilizer code. 

Lemma 10: If a pure ((n, K, R, d)) q Clifford subsystem 
code Q exists, then there exists a pure ((n, KR, d)) q stabi- 
lizer code. 

Proof: Let X be a classical additive subcode of F^" that 
defines Q, and let Y = X n X '. Furthermore, Theorem ^ 
implies that KR = q n /\Y\. Since Y C Y ±a , there exists 
an ((n, g n /|y|, d')) q stabilizer code with minimum distance 
d' = wt(y ±s — Y). The purity of Q implies that swt(y ±s - 
X) = swt(F- L =) = d. As Y C X, it follows that d! = 
swt(F ±s —Y) = swt(Y s ) = d; hence, there exists a pure 
((n, KR, d)) q stabilizer code. ■ 
As a consequence of the preceding lemma, it is straight- 
forward to obtain the following bounds on pure subsystem 
codes. 



Theorem 11: Any pure ((n, K, R, d)) q Clifford subsystem 
code satisfies KR < q n ~ 2d + 2 ^ 

Proof: By LemmaflOl there exists a pure ((n, KR, d)) q 
stabilizer code. By the quantum Singleton bound, we have 

KR < q n - 2d + 2 . m 

Corollary 12: A pure [[n, k, r, d]] q Clifford subsystem 
code satisfies k + r < n — 2d +2. 

Example 4 (Optimal Subsystem Codes): All the follow- 
ing codes constructed from Reed-Solomon codes over ¥ q i 
are pure and meet the bound in Theorem ^2 These codes 
are in that sense optimal subsystem codes. 



TABLE IV 
Optimal Pure Subsystem Codes 



Subsystem Codes 


Parent 
Code (RS Code) 


[[15,1,10,3H 4 
[[15,9,2,3]] 4 
[[15,10, 1,3]] 4 


[15,12,4J 4 2 
[15,4, 12] 4 2 
[15,3, 13] 4 2 


[[24,1,17,411s 
[[24,16,2,4]] 5 
[[24,17,1,4]] 5 
[[24, 19, 1, 3]]s 


[24,20,5] 5 2 
[24,5,20] 5 2 
[24,4,21] 52 
[24,3,22] 5 2 


[[48,1,37,6U7 


[48,42, 7| 7 2 



We can also show that the pure subsystem codes obey a 
quantum Hamming bound like the stabilizer codes. We skip 
the proof as it is along the same lines as Theorem II II 

Lemma 13: A pure ((n,K,R,d)) q Clifford subsystem 
code satisfies 

3=0 

V. SUBSYSTEM CODE CONSTRUCTIONS 

In this section, we give new constructions for pure subsys- 
tem codes. We begin with a proof of the simple, yet surpris- 
ing, observation that one can always exchange information 
qudits and gauge qudits in the case of pure subsystem codes. 

Lemma 14: If there exists a pure ((rt, K, R, d)) q Clifford 
subsystem code, then there also exists an ((n, R, K,>d)) q 
Clifford subsystem code that is pure to d. 

Proof: By Theorem [2] there exist classical codes 
D C C C W n q2 with the parameters {n,q n R/K) q2 and 
(n,q n / KR) q 2. Furthermore, since the subsystem code is 
pure, we have wt(D ±a \ C) = wt(D ±a ) = d. 

Let us interchange the roles of C and C ±a , that is, now 
we construct a subsystem code from C ±a . The parameters 
of the resulting subsystem code are given by 

((n, y/\D^\/\C^\, y/\C^\/\D\,wt(D^ \ C^))) q . 

We note that 

. V ggj \/\C^ | = tJ\C\IW\ = R and 

. y/\C^\/\D\ = y/\D^\f\C\=K. 
The minimum distance d! of the resulting code satisfies d! = 
wt(£'- La \C- La ) > wt(L»- La ) = d; the claim about the purity 
follows from the fact that wt(D ±a ) = d. ■ 

Before proving our next result, we need the following fact 
from linear algebra. 



Lemma 15: Let ¥ q be a finite field of characteristic p. Let 
C denote an additive subcode of F^™. There exists an F p - 
basis B generating the code C that is of the form 

B . . . , Z r , X r , ^r+lj • ■ • j Zr+j } 

where (x k \ x e ) s = = (z k \ z e ) s and (x k z t ) s = S k ,e- In 
particular, D = C (~l C " = (z r +i, ■ ■ ■ , z r+s ). It is possible 
to choose B such that it contains a vector z k of minimum 
weight swt(C). 

Proof: Choose a basis {z%, . . . , z r+ j} of a maximal 
isotropic subspace Co of C. If Co ^ C, then we can choose 
a codeword x\ in C that is orthogonal to all of the z k except 
one, say Z\ (renumbering if necessary). By multiplying with 
a scalar in F*, we may assume that (z\ | x 1 ) s = 1. If 
(Cq,Xi) ^ C, then one can repeat the process a finite 
number of times by choosing an x k that is orthogonal to 
{xi, . . . , Xk-i] until a basis of the desired form is found. ■ 

A subset {zk,Xk} of C with (z k \ Xk) s = 1 is called a 
hyperbolic pair. Thus, in the proof of the previous lemma, 
one chooses in each step a hyperbolic pair that is orthogonal 
to the previously chosen hyperbolic pairs. 

Theorem 16 ('Rain on your Parade Theorem'): Let ¥ q 
be a finite field of characteristic p. An ((n,K,R > l,d)) q 
Clifford subsystem code Q implies the existence of an 
((n,K,R/p,d)) q Clifford subsystem code Q s . If Q is 
exactly pure to d! , then the subsystem code Q s can be 
chosen such that it is exactly pure to d! as well. 

Proof: By Theorem \l\ there exists an additive code 
C < ¥ 2 q n with subcode D = C n C ±s such that K = 
q n /{\D\\C\) 1 l 2 ,R=(\C\/\D\) 1 / 2 ,d = swt{D^\C), and 
d! = swt(C). By Lemma IT31 one can find a F p -basis B of 
the form B = {zi,xi; . . . ; z r , x r ; z r+ i, . . . , z r+ j} such that 
(x k | xi) s = = (zk | z e ) s and (x k \ z e ) s = S k ,e- Notice 
that D = Cn C ±s = (zr+i, ■ ■ ■ , z r+j) by Lemma fT~5l 

Let C s be the additive subcode of C given by C s = 
span F (B \ {x r }). Then D s = C s H Cj-" = (z r , . . . , z r +j). 
It follows that |C S | = \C\/p and \D S \ = p\D\. Therefore, C s 
defines a subsystem code Q s = A S ®B S such that dim^4 s = 
q n /(\C s \\D s \) 1 / 2 = K and dimB s = {\C s \/\D s \f^ = 
Rip. 

Since D^r s C D ±s , any minimum weight codeword c £ 
Dj-' \ C a must be either in \ C or C. If it is in \ C, 
then swt(c) > d. If it is in C, then it is a linear combination 
of elements in B \ {x r }, since x r g" Djr 3 . This implies that 
c is contained in C s , contradicting our assumption that c is 
in D^ s \ C s . Therefore, swt(Djr 3 \ C s ) > d and we can 
conclude that Q s has minimum distance > d. 

For the purity statement, recall that D C D s C C s C C. 
The subsystem code Q is exactly pure to d! = swt(C). If 
swt(-D) = d! , then swt(C s ) = d!\ otherwise, swt(C \ D) = 
d' and we can choose z r+ i such that swt(z r +i) = d! . Then 
the subsystem code Q s is exactly pure to swt(C s ) = d! . ■ 

Corollary 17: An ((n, K, R, d)) q Clifford subsystem code 
that is exactly pure to d! implies the existence of an ((n, K, > 
d)) q stabilizer code that is (exactly) pure to d! . 

Proof: The corollary follows by repeatedly applying 
TheoremfT6lto the ((n, K 7 R, d)) q code and the derived code 



until the dimension of the gauge subsystem is reduced to one. 

■ 

We know that the MDS stabilizer codes arise from clas- 
sical MDS codes. In fact, the stabilizer code is MDS if 
and only if the associated classical code is MDS. We can 
therefore hope that good subsystem codes can be obtained 
from classical MDS codes. We show that the resulting 
subsystem codes must be pure. 

Lemma 18: If an ((n, K > 1, R > 1, d)) q subsystem code 
is constructed from an MDS code, then the resulting code is 
pure. 

Proof: Assume that C C F™ 2 is an [n, k,n — k + l] q 2 
code. If C ±a C C, then K = 1 contrary to our assumption. 
So assume that C ±a % C. Let k > n~k. Then D = CnC ±a 
must be smaller than C ±a . And dirn-D < n — k — 1. Hence 
wt(D ±a ) <(n-k-l) + l = n-k<n-k + l = wt(C). 
Hence the subsystem code is pure. Now assume that k < 
n-k. Now it is possible that C C C ±a . If C C C ±a , then 
R = 1. So C 1 C^. Now dim D < k — 1 from which it 
follows that wt(D J - a ) <k<n-k<n-k + l= wt(C). 
It follows that the subsystem code is pure. ■ 

VI. STABILIZER VERSUS SUBSYSTEM CODES 

In this section, we make a rigorous comparison between 
stabilizer codes and subsystem codes. Strictly speaking, 
subsystem codes contain the class of stabilizer codes; thus, 
in this section, we assume that the subsystem codes have a 
co-subsystem of dimension greater than 1. 

Clearly, there are difficulties in comparing the two classes 
of codes. Our "rain on your parade" theorem shows that 
Clifford subsystem codes cannot have higher distances than 
stabilizer codes. Their main edge lies in simpler error recov- 
ery schemes. We can quantify this in terms of the number of 
syndrome measurements required for error-correction. This 
is not necessarily the best method to compare the decoding 
complexity. However, it is certainly a reasonable measure 
if both codes use table lookup decoding. In the absence of 
any special algorithms for subsystem codes, we will proceed 
with this as the metric for comparison. 

A. Improving Upon Quantum MDS Codes 

In this subsection, we want to settle whether or not 
there exist subsystem code with parameters [[n,n — 2d + 
2, r > 0, d]] q . It turns out that the bounds that we have 
derived in Section IIVI will help in answering this ques- 
tion. Our best bounds are restricted to pure codes. Fortu- 
nately, it turns out that all subsystem codes with parameters 
({n, g"- 2 ^ 2 , R i d)) q are pure. 

Theorem 19: Any ({n,q n - 2d+2 , R,d)) q Clifford subsys- 
tem code is pure. 

Proof: If R = 1, then the claim follows from the fact 
that quantum MDS codes are pure, see [12]. 

Seeking a contradiction, we assume that there exists an im- 
pure subsystem code with parameters ((n, q"- 2d + 2 ; d)) q , 
exactly pure to d! < d and R > 1. It follows from Corol- 
lary [n] that it is possible to construct a stabilizer code with 
distance > d that is (exactly) pure to d! . Then the resulting 



stabilizer code has the parameters ((n,q n ~ 2d+2 ,d)) q and 
is impure. But we know that all quantum MDS codes are 
pure [12], see also [5, Corollary 60]. This implies that 
d! > d contradicting the fact that d' < d\ hence, every 
((n, q n ~ 2d + 2 ^ R y d)) q code is pure. ■ 

The next theorem explains why Poulin did not have any 
luck in finding an [[5, 1, r > 0, 3]]2 subsystem code. 

Theorem 20: There do not exist any Clifford subsystem 
codes with parameters ((n,q n ~ 2d+2 ,R > l,d)) q . In partic- 
ular, there do not exist any [[n,n—2d+2, r > 0, d]]„ Clifford 
subsystem codes. 

Proof: Seeking a contradiction, we assume that a 
subsystem code with parameters ((n, q n ~ 2d+2 , R > l,d)) q 
exists. By Theorem 1191 an ((n 7 q n ~ 2d+2 ,R,d)) q subsystem 
code must be pure. It follows from Theorem that a pure 
subsystem code with these parameters must satisfy 

n-2d+2ji < _n-2d+2 

Therefore, we must have R = 1, contradicting our assump- 
tion R > 1. ■ 

B. Better Than Quantum MDS Codes 

In this subsection, we compare once again quantum MDS 
stabilizer codes against subsystem codes. We require that 
both codes are able to encode the same amount of infor- 
mation and have the same distance. However, this time, 
we do not restrict the length of the codes. Our goal is to 
determine whether the subsystem code can improve upon an 
optimal quantum MDS stabilizer code by fewer syndrome 
measurements. 

We insist that the codes are F 9 -linear, since in this case the 
number of syndrome measurements can be directly obtained 
from the code parameters. Indeed, recall that an F g -linear 
[[n, k, r, d]) q subsystem code requires n — k — r syndrome 
measurements, and an F g -linear \[n' , kl , d'\\ q stabilizer code 
requires n' — k' syndrome measurements. 

Theorem 21: If there exists an F g -linear [[k+2d— 2, k, d]] q 
quantum MDS stabilizer code, then an F g -linear [[n, k, r, d]] q 
subsystem code satisfying 

k + r<n-2d+2 (1) 

cannot require fewer syndrome measurements than the sta- 
bilizer code. 

[We remark that any pure [[n, k, r, d]] q subsystem code 
satisfies the inequality Q by Theorem II 111 

Proof: Seeking a contradiction, we assume that the 
subsystem code requires fewer syndrome measurements than 
the quantum MDS code, that is, we assume that k + 2d — 
2 — k > n — k — r. This implies that k + r > n — 2d + 2, 
contradicting our assumption that k + r < n — 2d + 2. ■ 

Now, we can partially answer the question when an ¥ q - 
linear [[n, k, r, d]] q subsystem code will lead to better error 
recovery schemes than the quantum MDS codes. 

Corollary 22: Suppose that an F g -linear \[k + 2d — 
2, fc, d]] q quantum MDS code Q exists. Then an F g -linear 
[[n, fc, r, d]] q subsystem code that beats the the stabilizer code 
Q must be impure and must satisfy k + r > n — 2d + 2. 



Proof: We know from Theorem that all pure 
[[n, k 1 r, d]] q codes satisfy k + r < n — 2d + 2. But Theo- 
rem implies that such a code cannot have fewer syndrome 
measurements than the F g -linear MDS code. Hence, the 
subsystem code, if it is better, must be impure and have 
k + r > n - 2d + 2. ■ 

C. Better Than Optimal non-MDS Stabilizer Codes 

We know that MDS codes do not exist for all lengths, 
so it is reasonable to consider optimal stabilizer codes that 
are non-MDS. In this case, the comparison is slightly more 
complicated. An [[n, k, r, d]] q subsystem code could be better 
than an optimal [[n 1 , k, d]} q stabilizer code. That in itself does 
not guarantee that the class of subsystem codes is superior 
to the class of stabilizer codes. 

For instance, the shortest code to encode 2 qubits with 
distance 3 is [[8, 2, 3]] 2 (see [3]). Suppose that an [[8, 2, 1, 3]] 2 
code exists. This subsystem code requires only 8 — 2 — 
1 = 5 syndrome measurements as against the 8 — 2 = 6 
measurements of the optimal stabilizer code. To conclude 
that the subsystem codes are better than stabilizer codes 
would be premature, for there exists an [[8, 3, 3]] 2 code 
(cf. [3]) that requires 8 — 3 = 5 syndrome measurements 
and encodes one more qubit than the subsystem code. It is 
therefore necessary to compare the subsystem code with all 
optimal [[n',k' > k,d]] q stabilizer codes, where n' ranges 
from n — r to n. Only if the subsystem code requires 
fewer syndrome measurements in each case, then we can 
conclude that the class of subsystem codes leads to better 
error recovery schemes. 

We do not know precisely the properties of such subsystem 
codes. For instance, we do not know if such subsystem code 
is required to be impure or if it must satisfy k+r > n—2d+2. 

Next, we turn our attention to a slightly different question 
which shows that in general whenever good subsystem codes 
exist, good stabilizer codes also exist. 

D. Subsystem Codes and Stabilizer Codes of Comparable 
Performance 

The reader will perhaps wonder why one cannot simply 
discard the gauge subsystem to obtain a shorter quantum 
code without sacrificing distance or dimension. The reason 
why we cannot do so is because, in general, there is no 
one to one correspondence between the gauge qudits and the 
physical qudits. Yet, our intuition is not entirely misguided 
as the following result will show. 

Theorem 23: Let ¥ q be finite field of characteristic p and 
1 < q k < q n a power of p. Let r be an integer such that 
< r < n, and 

d-i , , 

£ (<? 2 - iy(q n+k+r - q n+r - k ) < (p - l)(q 2n - 1) 
3=1 ^ 

holds, then there exist both an ((n,q k ,q r ,> d)) q Clifford 
subsystem code and an ((n — r, q k , > d)) q stabilizer code. 



Proof: By hypothesis 

E0;y-i)'<(p-i) 



(q 2n - 1) 



holds and Theorem [7] implies the existence of an 
((n, q k ,q r ,> d)) q Clifford subsystem code. We can rewrite 
the RHS of the inequality as 



RHS 



(P-1) 
(P-1) 
(P-1) 





- 9" 


- 71 — V 


q k 


- 9" 


-k 


q n-r 


- 9" 


-n-\-r 
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- <T 


-k 


q n-r 
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-n+r 


qk 


- 9" 


-k 



7 — n+r 



+ (p - 1) 



— n+r 



i-k 



i-k 



<1, if r>0 



Now under the assumption r < n, we obtain a lower bound 
for LHS as follows. 



d-l 



;> 2 -i)+e 

d-l 

= E 



d-l 

^E 



(? 2 - 1) 3 ' 



( g 2 -ir"+r( g 2 -l), 



(g 2 - 1)J = LHS. 



Since we know that LHS < RHS we can write 

d-l 



(g 2 - iy + r(q 2 - 1) 




d-l 

E 



\ q 2 -iy < p-i - 

J 



q k — q k 



Then by Corollary |S] there exists an ((n — r, q , > d)) q 
stabilizer code. ■ 
While they might differ in their distance, the preceding 
theorem indicates that in many cases, whenever a good 
subsystem code exists, then there will also exist a good 
stabilizer code encoding as much information and having 
comparable distance and of shorter length. The assumption 
of integral r may not be much of a restriction in light of 
Theorem [16] 



In comparing the complexity of the error recovery schemes 
for the two codes, we run into a small problem since we 
do not know if the codes are F g -linear. Actually, if we use 
the stronger result of [5, Lemma 31] and insist that n = 
k mod 2, then we can show that the stabilizer code is ¥ q - 
linear. This guarantees that the stabilizer code will require 
n—k—r syndrome measurements which is comparable to that 
of an F 9 -linear subsystem code. It appears then, quite often, 
subsystem codes do not offer any gains in error recovery, as 
there will exist a corresponding stabilizer code that encodes 
as many qudits, of similar distance and equal complexity of 
decoding. 

VII. CONCLUSION 

In this paper we investigated subsystem codes and their 
connections to classical codes. We derived lower and upper 
bounds on the parameters of the subsystem codes. We settled 
the question whether or not there exist n — 2d + 2, r > 
0,d]] q subsystem codes exist. We showed that pure F g - 
linear subsystem codes do not lead to any reduction in 
complexity of error recovery as compared with an F 9 -linear 
MDS stabilizer code of equal capability. As a consequence 
we concluded that the subsystem codes that outperform the 
quantum MDS codes must be impure. Further, we showed 
that quite often the existence of a subsystem code implies 
the existence of a stabilizer code of comparable performance 
and complexity of error recovery. 
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